Ground state non- universality in the random-field Ising model 
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Two attractive and often used ideas, namely universality and the concept of a zero temperature 
fixed point, are violated in the infinite-range random- field Ising model. In the ground state we show 
that the exponents can depend continuously on the disorder and so are non-universal. However, we 
also show that at finite temperature the thermal order-parameter exponent 1/2 is restored so that 
temperature is a relevant variable. The broader implications of these results are discussed. 
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Co-operative behavior in disordered systems can usu- 
ally be concisely characterised using scaling theories. 
These scaling theories contain scaling exponents and it is 
particularly important and satisfying if these exponents 
are independent of the fine details of the model, that is 
they are in some sense "universal" . If this holds it allows 
the theorist to study the simplest or most convenient 
model in a class in order to find the scaling exponents 
and, more importantly, that experiments should show 
the same exponents as the theory even though they may 
look very different on short length scales. Universality 
has been spectacularly successful in the study of phase 
transitions as a function of temperature, culminating in 
the development of the renormalisation group jjj . There 
the critical exponents, usually, only depend on the sym- 
metry of the order parameter and the spatial dimension. 
Due to the fact that scaling theories also work in many 
disordered systems it is natural to try to extend the ideas 
which work so well for thermal phase transitions to the 
disorder case. 

Universality in disordered systems assumes that the 
critical exponents should not depend of the type of dis- 
order, provided the disorder distribution is short-range 
correlated and provided it is not too broad. Universality 
with respect to disorder has been confirmed in some sys- 
tems, with notable examples being percolation 0] and the 
problem of a directed polymer in a random medium 
However universality in disordered systems is unproven 
in general and indeed it has been questioned in the spin 
glass problem where there appears to be qualitative dif- 
ference between the behavior in the presence of Gaussian 
as compared to bimodal disorder Q] . More recently uni- 
versality has even been questioned in the random field 
Ising model ||[| , which is one of the simplest models of 
a disordered material. Moreover the experimental tests 
of the random field Ising exponents rely on universality 
as the experiments are carried out on diluted anti- 
ferromagnets in a field Q which are expected to lie in 
the same universality class (these experiments are also 



plagued by kinetic effects due to the large barriers which 
exist in random magnets). We show that universality 
fails in the ground state of the infinite-range random- field 
Ising model, as the exponents may vary continuously with 
the type of disorder. 

However, we also show that universality is restored at 
any finite temperature in the sense that at finite temper- 
ature the order parameter exponent is always 1/2, when 
the transition is continuous. This implies that another 
important concept in disordered systems, the concept of 
a zero temperature fixed point, is violated in this model. 
The origin of the concept of a zero temperature fixed point 
is that disorder usually provides a stronger perturbation 
than thermal fluctuations (which in turn are usually a 
stronger perturbation than quantum fluctuations) . Thus 
a study of the ground states of disordered systems can 
lead to scaling theories which are qualitatively correct at 
finite temperatures. This is particularly attractive since 
there now exist methods for finding the exact ground 
states of many quenched random systems JToj ] . In the 
vernacular of random systems it is often stated that there 
exists a zero-temperature fixed point which controls the 
behavior at finite temperature. In particular, in both the 
random field Ising model 1 1 -|l4|] and in spin glasses [|J, 
scaling theories are frequently based on the assumption of 
a zero-temperature fixed point. However, we show that 
the mean-field theory of the random field Ising model 
is not, in general, controlled by a zero-temperature fixed 
point. 

We first demonstrate that the critical exponents can 
take on a range of values in the ground state of the 
random-field Ising model. The Hamiltonian for this 
model is, 



H = - Jo S i S i ~ Yl hiSi = NE * 



NE 



f 



(1) 



where the first sum is over all spin pairs and Jo = J/N 
where N is the number of sites in the lattice to ensure 
an extensive energy. When the distribution of random 
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fields is narrow, the exchange term dominates and the 
system is a ferromagnet (in dimensions greater than or 
equal to three), while when the random field distribution 
is broad the random field dominates and the system be- 
comes a paramagnet. We take distributions of random 
fields which have mean zero and width 5h, and consider 
the key ratio H = 5h/J which measures the strength of 
the random field in comparison to the exchange. 

Consider a spin subspace in which the magnetisation, 
to, is fixed, ie. to = (n + — nJ)/N, where n + is the num- 
ber of up spins in the configuration and n_ is the number 
of down spins. It is easy to find the lowest energy state 
for fixed to. The exchange energy is given by, 



Eex (m) 



-J 



2N 2 



2n + n_) 



-Jto 



(2) 



Due to the fact that the exchange is of infinite range, 
all configurations at fixed to have the same energy and 
so are combinatorially degenerate. The field term splits 
this degeneracy by choosing the configuration which has 
the smallest field energy. This is achieved by satisfying 
the largest random fields and leaving the smallest possi- 
ble fields unsatisfied. If the distribution of random fields 
is P(h) (which we assume to be symmetric about the 
origin), then in the large lattice limit, we have, 



E f = -2 



<j] = -i 1 dhhP(h) + 

lo Jo 



h c (m) 



dhhP(h) (3) 



where the first term is the ideal field energy in which ev- 
ery spin is oriented in the direction of its local field, and 
the second term is the energy cost due to the fraction of 
fields which are unsatisfied. The fraction of fields which 
are unsatisfied is determined by the magnetisation, 



to = 2 



h c (m) 



P(h)dh. 



(4) 



The ground state is found by determining the value of to 
which minimizes the energy (1-3), given the constraint 
(4). Carrying out the variation yields, 



° {E - + E A = - Jto + 2h c (m)P(h c (m)) dkc{m) 



dm 



dm 



(5) 



By taking a derivative of Eq. (4) with respect to to (us- 
ing the chain rule) we find, 1 = 2P(h c (m))dh c (m) / dm 
Using this to remove dh c (m) / dm from the RHS of Eq. 
(5) and setting Eq. (5) to zero, we find that the cutoff 
field is related to the magnetisation via h c (m) = Jto. 
Subsitution of this into Eq. (4) yields the ground-state 
mean-field equation 



2 I P(h)dh. 
o 



(6) 



This equation gives the magnetisation values at which 
the energy is extremal. Note that to = (the param- 
agnet) is always an extremum, as expected. Since we 



are treating the case of symmetric random fields, we can 
restrict attention to the case where < to < 1 . To deter- 
mine whether an extrema is a maximum or a minimum, 
we need to evaluate the curvature near the extremum, 



d 2 (E ex +E f ) 
dm 2 



= — J 



1 



2P{h c {m s )) 



(7) 



Finally, in order to determine the ground state, we need 
to compare the free energies of the solutions to the mean- 
field equation (4) with the energy of the magnetised state 
ie. to = 1. 

An elegant result due to Aharony |l9] states that the 
nature of the finite temperature phase transition in the 
infinite range random field Ising model depends on the 
curvature of the disorder distribution at the origin. Bi- 
modal distributions lead to a first order jump in the 
order parameter at low temperatures (and hence a tri- 
critical point at finite temperature), while unimodal dis- 
tributions exhibit continuous transitions with exponent 
(3 = 1/2, as originally found by Schneider and Pyttc llfl] 
for the case of Gaussian disorder. However, we now show 
that the exponent (3 = 1/2 is not universal in the ground 
state. 

We show that may change continuously with the 
disorder, by considering the distribution of random fields 
given by, 

p(h) = - (\hl)v) -H<h<H 1 (8) 

with y > 0. In the limit y — > oo P(h) — > uniform so that 
a first order behavior is expected, while if j/ > 2 it looks 
like a Gaussian near the origin so we expect a continuous 
transition with j3 = 1/2. In the following discussion, we 
take J = 1, so that H has been normalised by J. The 
distribution (8) is the first two terms in the expansion 
of the stretched exponential, exp(— (\h\/H) v ) which has 
the same critical behavior as (8). However, Eq. (6) can 
be solved exactly for the case (8) to yield (in addition to 
to = 0), 



H(y(H c - H)) 1/y for 1 < H < H c 



(9) 



where the critical field is given by, H c = (y + l)/y. The 
lower bound on H is due to the cutoff in Eq. (8). For 
H < 1, the exchange always wins and the magnetisation 
is to = 1. From (7), the second derivative is, 



d 2 (E ex +E f ) 
dm 2 



= -1 



yH 



(10) 



Thus the curvature at zero magnetisation changes from 
positive (a minimum) for H > H c to negative for H < 
H c . It is also easy to show that the solution m s is always 
a minimum. Evaluating the energies at the three solu- 
tions to = 0, to — m si m = 1 yields a behavior typified 
by Fig. la. For H > H c , the ground state has to = 
and the system is a paramagnet, for 1 < H < H c , the 
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ground state has m = m s and is magnetised, while for 
H < 1, the magnetisation saturates. This behavior is 
summarised in the phase diagram of Fig. lb. The criti- 
cal exponent (3 = 1/y on this upper curve in this figure 
and is clearly non-universal in the ground state. 

However when the transition is continuous, Aharony 
[jl5j has demonstrated that at finite temperature (3=1/2 
based on the mean-field equation for the random-field 
Ising model, 



dhP(h)taah(m/T + h/T). 



(11) 



We now reconcile the ground state result (9) found above 
with the finite temperature behavior found from equa- 
tion (11). Assuming that P(h) is symmetric, equation 
(11) can be reduced to, 



2tanh(2m/T) 



P{h)dh 



1 



cosh(2fe/T) \ ' 
cosh(2m/T) / 



(12) 



From this expression, it is seen that there are two 
regimes, m/T >> 1 and m/T << 1. At zero temper- 
ature only the first regime holds, while at any finite tem- 
perature the second regime is applicable very close to the 
critical point. 

When m/T » 1, tanh(m/T) -> 1 and 
cosh(2ft,/T)/cosh(2m/T) -> exp(2(/i - m)/T), which 
yields, 



m = 2 



P(h)dh 



1 + exp(2(/i - m)/T) 



m/T — > oo. (13) 



Now note that this expression looks like a Sommcrfeld in- 
tegral for the free fermi gas, with the Fermi energy given 
by, ef — m. The leading term at low temperatures is 
then the integral of P(h) up to the Fermi energy, and 
hence is equivalent to the ground state result given in 
Eq. (6). 

However at any finite temperature, there is a regime in 
which the magnetisation is small compared to the tem- 
perature, m/T « 1. In that case, cosh(2m/T) — > 1, 
and equation (11) reduces to the mean field theory for 
the thermal transition, but with a renormalised coeffi- 
cient which depends on the field distribution, ie., 



2I(H, T)tanh(2m/T) m/T -> 



where, 



I(H,T) 



P(h)dh 



l + cosh(2h/T)' 



(14) 



(15) 



Note that there is a factor of two difference in the ar- 
gument of the tanh as compared to the thermal mean 
field theory. However the critical temperature and crit- 
ical exponent are the same. For any finite temperature, 
provided m << T, an expansion to third order in m of 
Eq. (14) shows that the magnetisation approaches zero 



with exponent (3 = 1/2. Moreover, Eq. (14) with (15) 
shows that the critical field and temperature are related 
to each other through the relation, 



T = iI(H c {T),T), 



(16) 



provided the magnetisation is continuous at the transi- 
tion. For the probability distribution given in Eq. (8), 
this reduces to, 



Hc(T) u HcP) 



T 

— ) y 



x v dx 



1 + cosh(x) 
(17) 



The magnetisation as a function of field is given in Fig. 
2a for y = 1/2. From this figure it is seen that the critical 
exponent in the ground state is different than that at fi- 
nite temperatures, and clearly illustrates the fact that 
temperature is a relevant variable. The temperature- 
field phase diagram is presented in Fig. 2b for the two 
cases y = 1/2 and y = 2. There is a sharp shift in the 
phase boundary with temperature for cases where y is 
small (rapidly decaying field distributions near the ori- 
gin), which is strong indicator that temperature is rele- 
vant. 

We have demonstrated the failure of universality in the 
ground state of the mean field theory of the random field 
Ising model. In addition the concept of a zero temper- 
ature fixed point is invalid. The fact that the Gaussian 
distribution of random fields does have exponent (3=1/2 
in the ground state is atypical and should not be ex- 
pected unless the disorder distribution is quadratic near 
the origin. Finite temperature introduces thermal fluctu- 
ations which are also Gaussian which is the reason that 
the Gaussian distribution of disorder is special and atyp- 
ical. 

At first blush, our results raise serious questions about 
scaling theories of disordered systems based on a zero 
temperature fixed point, and about the applicability of 
numerical studies in the ground state to finite temper- 
ature properties. However this may not be the correct 
conclusion. Instead the conventional mean-field theory 
described here may be pathological and not typical of 
the behavior in finite dimensions. That itself would be 
a rather surprising result, which for example could be 
due to a renormalisation of the disorder distribution to 
a Gaussian under rescaling in finite dimensions. These 
issues can only be resolved by careful studies of univer- 
sality to disorder in finite dimensions, which is a difficult 
task except at zero temperature where exact numerical 
calculations are possible. 

This work has been supported by the DOE under con- 
tract DE-FG02-90ER45418. 
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FIG. 1. (a) The ground-state energy as a function of the 
width of the magnetic field distribution (H), for the case 
y — 1/2. The flat curve is for m — 1, the linearly decreasing 
curve is for m = while the third curve is for the solution 
m 3 given in Eq. (9) of the text, (b) The upper curve is the 
dependence of the critical field on the exponent in the field 
distribution (i.e. H — (y + l)/y). Above this line the mag- 
netisation is zero. The lower line is H = 1, below which the 
magnetisation is saturated (i.e. m — 1 for H < 1). Between 
these two lines the magnetisation obeys Eq. (9) of the text, 
with the magnetisation going to zero with exponent 1/y at 
the upper curve. 
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FIG. 2. (a) The magnetisation as a function of the width 
of the random field (H) for y = 1/2, for two temperatures. 
The upper curve is for zero temperature, while the lower 
curve is for T = 0.2. (b) The H — T phase diagram (T vs. 
H c (T)/H c (0)) for y = 2 (upper curve) and y = 1/2 (lower 
curve) found from solving Eq. (17) of the text. 
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